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We obtain effective potential of 0(A^)-symmetric (f)^ theory for large starting with a finite lattice 
system and taking the thermodynamic limit with great care. In the thermodynamic limit, it is 
globally real-valued and convex in both the symmetric and the broken phases. In particular, it has 
a fiat bottom in the broken phase. Taking the continuum limit, we discuss renormalization effects 
to the flat bottom and exhibit the effective potential of the continuum theory in three and four 
dimensions. On the other hand the effective potential is nonconvex in a finite lattice system. Our 
numerical study shows that the barrier height of the effective potential flattens as a linear size of the 
system becomes large. It decreases obeying power law and the exponent is about —2. The result is 
clearly understood from dominance of configurations with slowly-rotating field in one direction. 
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1 Introduction 



Effective potential is widely employed in various contexts of physics for studying phase transition. 
In the language of statistical physics, it is a free-energy density with some dynamical variable fixed. 
The fixed dynamical variable is no longer "dynamical", and is called an order parameter. A value 
of the order parameter in the thermodynamic limit specifies a phase of a model. As shown in later 
section, the value of the order parameter is dynamically determined in a thermal equilibrium such 
that it gives a global minimum of the effective potential. Therefore we can recognize which phase is 
realized by looking at a form of the effective potential. 

For example, let us consider a classical ferromagnetic spin system. The order parameter in this 
case is a component of the spin variable averaged over all lattice sites. According to the Landau- 
Ginzburg theory of phase transition, the effective potential is phenomenologically given as a poly- 
nomial of the order parameter: it has the unique minimum at the origin in the high-temperature 
(symmetric) phase. That is, a value of the order parameter is zero in this phase. As temperature 
is lowered, the effective potential continuously changes and it eventually comes to have form of a 
double-well. Namely the order parameter does not vanish any more and we recognize that the system 
is in the low-temperature (broken) phase. 

Although the Landau- Ginzburg theory gives a good qualitative picture of phase transition, the 
form of the phenomenological effective potential is not globally correct because the effective potential 
of the spin system must be convex even in the low-temperature phase. More precisely, as proved 
by O'Raifeartaigh et al., the effective potential of a lattice scalar model with a nearest neighbor 
interaction plus on-site potential must be always convex even if the classical potential is not[jll. It 
means that a form like a double- well is ruled out. 

How does the effective potential look like in the low-temperature phase? Numerical simulations for 

theory on a lattice suggest that it has a flat bottom |l| . Similar effective potential is analytically 
obtained by LangerjQ in the spherical model with a help of the so-called sticking argument [pO| , 
The fiat region reflects dominance of configuration with domain structures [^, or dominance of 
spatially slowly varying configuration |ll| . 

A usual way of computing the effective potential of (p'^ theory is a loop expansion from a uniform 
background, which is considered to be valid in a weak coupling region[||, It gives the classical 
potential in the leading order. However, when the classical potential is not convex, the loop expansion 
develops an imaginary part due to unstable modes of fiuctuation0. Systematic formulations that 
circumvent the unstable fluctuation have been proposed and one can obtain the (almost) convex 
effective potential by these methods|P, |ll], [T^. However, these results do not give us quantitative 
understanding because one cannot explicitly find any controllable expansion parameters. Necessity 
of non-perturbative and quantitative approach naturally leads us to investigating the case where the 
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number of fields, N, is sufficiently large and to carrying out the 1/N expansion. 

Investigation of phase structure of the large- 0^ theory using the effective potential has been 
extensively carried out for more than two decades |T^, |T^, ^, |TB|, 17]. In order to determine a 
value of the order parameter by the effective potential, we need to know where a global minimum 
is. However, most of the authors encountered trouble, i.e., an imaginary part appears in the broken 
phase and they fails to obtain the effective potential globally. In our knowledge, only exception 
is ref.0, where a non-uniform external field is used as a "regularization" . Although the effective 
potential seems to be globally obtained in that paper, the ansatz used to solve a saddle-point equation 
is incorrect. This point will be discussed more detail in sect. 5. 

Thus, we attempt to find a global form of the effective potential of the 0^ model in the large- 
limit, which is the main subject of this paper. 

We solve the saddle-point equation with giving a special attention to asymptotic property of 
the solution when the volume of the system tends to infinity. If one solves it after taking the 
thermodynamic limit, there are no real solutions for some range of order parameters in the broken 
phase. This is the reason why the appearance of the imaginary part explained above. On the other 
hand, if one solves it in a finite system and then takes the thermodynamic limit, one can obtain a 
globally real-valued effective potential. The idea of going back to a finite system when one solves a 



saddle-point equation was used to reduce mathematical manipulation of the sticking argument |]TB|, |K 

The rest of the paper is organized as follows. In sect. 2, we define 0(A^)-symmetric 0^ model 
on a lattice and give a formulation of the 1/A^ expansion of the effective potential in a finite lattice 
system. In sect. 3, we study asymptotic property of a solution of the saddle point equation derived 
in sect. 2 when the number of lattice sites of the system tends to infinity. Then we exhibit forms 
of the effective potential in the thermodynamic limit for various dimensions. Also, we numerically 
solve the saddle-point equation in a finite system and quantitatively study how effective potential in 
a finite system approaches to that in the thermodynamic limit. In sect. 4, we consider the continuum 
limit of the model and discuss the renormalization effects to the effective potential in three and four 
dimensions. Section 5 is devoted to summary and discussion. 

Here we would like to comment on definition of effective potential. The conventional way to 
define the effective potential is to introduce an external source coupled to some dynamical variable. 
The expectation value of the dynamical variable, which is the order parameter in this case, and also 
the free-energy density, are a function of the external source. We can regard the free-energy density 
as a function of the order parameter by the Legendre transformation]^. The free-energy density 
constructed in this way is called the effective potential. 

Alternatively we can fix a dynamical variable by imposing a constraint on a measure of the 
partition function so that the dynamical variable takes an given value (see eq.(|D). The free-energy 
density derived from the partition function using this measure is also called the effective potential 
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(or called the constraint effective potential |]T|). 

The relationship between the two is clarified in refs.[|5|, |l|: the former is the convex hull of the 
latter. O'Raifeartaigh et al. also show that the effective potential of the lattice (f)'^ theory derived 
from the second definition is convex. Therefore the two definitions give the same effective potential 
in the case of the lattice (p^ theory. In this paper we adopt the second definition. 



2 The expansion 

We study the 0(A^)-symmetric 0^ model on a D- dimensional hypercubic lattice. The lattice spacing 
is set to be unity. There are A^-scalar fields 0", a = 1, ■ ■ ■ , N defined on each lattice site x. We impose 
periodic boundary conditions: 0" = (p^+Lei^: = ■ ■ 'i-D, where e'^ is the unit vector pointing to 
the /i-direction and L is an integer standing for the length of the side of the hypercubic lattice. The 
model is defined by the action 

So = Eil^i-^ + + ^{€r.f}^ (1) 

where the summation over a is implicitly taken. The lattice Laplace operator A acts on fields as 

A0: = E(0S+e, + 0^e,-20S)- (2) 

Since we are interested in the broken phase, ttIq is set to be negative, while qq is taken to be positive 
because of the requirement that Sq should be bounded from below. 



We can replace Sq with the following action by introducing an auxihary field Xx as follows jlS 



1 N Nm"^ 

S = E{^0S(-A + XxWx - ^Xl + -^Xx}. (3) 
X ^ ^9o 9o 

In fact, the integration over Xx in the partition function is easily performed and the resulting action 
equals 5*0 up to a constant. Note that the sign of the xl term is negative: Xx is integrated along the 
imaginary axis. 

The effective potential U of the finite system is defined byp2|, |T3| || 

exp{-nuin,cp')}^ [Hdrxm-^^n UdxxeM-s}, (4) 

x,a X 

where (p"^ = (p°'ip°' and Q is the total number of the lattice sites. We choose 

a = l,---,iV (5) 
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as order parameters. It is obvious that f/ is a function of due to the 0{N) symmetry. 
The partition function Z of the system is written in terms of U{Q, tp^) as 

Z = I []rf(/^° exp{-fi U{n,ip^)]. (6) 

In the limit of — oo, the integration on the right-hand-side of eq.(^ is evaluated by a global 
minimum of U. i.e., the global minimum of f/(oo,<^^) determines a value of the order parameter 
realized in the thermodynamic limit. If there are multiple global minima, we need to choose one of 
them for looking at spontaneously symmetry breaking. A practical way of doing this is presented in 
sect. 3 (see also [|l|). A physical meaning of this procedure is explained in the literature [0, p^] , 
for example. 



\<4>l<t>l + ^{€€? (7) 



The classical potential 



takes the minimum on the (A^ — l)-dimensional sphere with the radius \J —2NmQ/ qq in </)-space(the 
target space). We expect that the values of Nip^ minimizing U lies in a range Nif^ < 0{N) as well as 
those of the classical potential, so that we shall search a minimum of U in that range: ip^ < 0{N^). 

The main purpose of this section is to derive a formulation of the 1/A^ expansion of U{Q,ifi^) 
following the spirit of refs.[^ That is, we perform integration over (p^ and obtain a theory 

described by Xx- The integration over Xx is approximated by the saddle point method, which is 
applicable for large A^. 

In order to carry out this program, we use an integral representation of the ^-function in eq.(^ 

n m - v^^-^) = £2 2^ ^^p{^'^(E € - nVNvi}- (8) 

We can readily integrate over 0". and r]"-. The result is, up to a constant, 
exp{-Q U{Q,ip^)}= [l[dXxexp{-NS,}, 

X 

1 1 1 TO^ 1 

[x] = -Tr ln(-A + x) + ^ In / + E(-^X^ + —Xx) + ^^T'^fi'- (9) 



Here we defined the operator x 
and the functional / 



xV) = Xx\x), (10) 



x^y '-^ ~r A 

The next step is to find the solution of the saddle-point equation 

\y\ , , 

— = 0. 12 

^Xx 
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We assume translationally invariant solution x* ■ In this case, the above equation (|T2|) is reduced to 

The solution x* is regarded as a function of ip"^ and Q. We get the leading-order effective potential, 
Uo{Q, ip^), by inserting x* into Si[x\: 

= :^Eln(-A(p) + x*)-— + '-X* + ^^'X*- 14 

2iip_^o 25-0 fi-o 2 

Further investigation to the leading order will be extensively carried out in later section. 

In the following we formulate calculation of higher-order corrections. The 1/N expansion is the 
expansion around the saddle point. Writing 

Xx = X* + ^^^^^^ (1^) 

then we expand 5*1 in terms of ^x'- 

NS,[x] = NS^lx*] + - T^TT^lxlUy + E ry tt 'j^lxlU ■ ■ ■ W 

Note that a linear term in C,x vanishes due to the saddle-point equation. 

The explicit form of the quadratic term is computed using the following formulae: 

Trln(-A + x) = E M-A(p) + X*) + ^ E f^o(p) E + ^ E ex4n(a; - 2/) + ■ ■ ■ , 

p ViViZ p X x,y 

Here A{p) is a Fourier component of the Laplace operator given by 

D 

A(p) = 2(Ecosp^-Z}), (18) 
11=1 



and 



Goix) = ^Ee^'^^o(p), Goip)^ , 
p -A{p) + x* 

n(^) = ^Ee^'^'n(p), m^^j:Go{p-q)Go{q). (19) 
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The result is 

D-\P) = ^(n(p)-^^Po) + ^ + (<^^-^)(Go(j>)-^5po). (20) 
We can read the propagator of from the quadratic term: 

{Uy) = ^T.^'''^^''^D{p). (21) 

Integrating over C,x with ignoring the cubic and higher-order vertices in (|T6|) gives the next-to- 
leading correction, which is of order N^. Higher-order corrections are calculated in a diagrammatic 
way. We see that the p-point vertex (p > 3) is suppressed by the factor of A^~p/2+^. Therefore a loop 
diagram that consists in pi-, ■ ■ -, p^- point vertices contributes to the order of A^~(pi+ "+Pfc)/2+'=. Xhe 
first nontrivial two-loop diagrams are generated by one 4-point vertex, or two 3-point vertices, which 
gives correction of order of A^~^. 

Finally we go back to the integration over (/)° that leads to eq.(^ and comment on an allowed 
range of x*- The first term of Uq{Q, y?^) in eq.(0) is derived from the integration over the oscillation 
modes of 0". The condition that the integration should converge is that 

5 + Re X* > 0, (22) 

where S is the smallest eigenvalue of —A in the oscillating modes: 

5 = 2(l-cos(— )). (23) 

3 Solution of the saddle point equation and a form of the 
effective potential 

3.1 Effective potential in the thermodynamic Umit 

In this subsection we shall study the leading-order effective potential of the lattice 0'* theory in 
the thermodynamic limit Q ^ oo. In particular, we would like to clarify the form of the effective 
potential in the broken phase. We denote 

u,in,^')^^Uoin,<p'), (24) 



6 




y = Hx) 



y = f + 2ml/ go 



Fig. 1. The graph of h{x) in £> = 4, where n = 12^. The 
couphng constants are go ~ 80, ^ —15. In this case, 
—S « —0.26. There exists unique solution x* for all t. 



for convenience. Since Mo(f2, 9?^) is an 0(A^)-symmetric function, we can put if^ = t, ip^ = {2 < 
k < N) without loss of generality. From eqs. (]TB|) and we have 

u',{n,f)=X% (25) 

u'^{Q,t') = x* + X*'t, (26) 

where ' means differentiation with respect to t. We wish to know asymptotic behavior of x* when 
Q — > 00. Let 

h{x) = -^y] + -X. (27) 

n^,-Aip) + x go 

Namely, h{x*) is the right-hand-side of the saddle-point equation(^). The solution x* is read from 
the crossing point having the graph of y = h{x) and a horizontal line representing y = + 2m1/ g^, 
as is shown in Fig. 1. 

Since h{x) is strictly increasing and takes any real value for x G (—5, 00), there is a unique real 
solution satisfying (|22|) for arbitrary t and Vt. The sign of x* depends on whether the horizontal line 
lies above or below /i(0). If x* is positive, the summation over p in h{x*) becomes the integration in 
the thermodynamic limit: 

lim hix*) = - I '^''^^ . ^ + -x\ (28) 

^ J (27r)^-A(p) + x* % ^ ' 

where the momentum {1 < ^ < D) runs from — tt to tt. The saddle-point equation ([l3|) in this 
case was essentially studied in ref. [0. On the other hand, if x* is negative, we cannot employ 
eq.(p8D any more and have to take the thermodynamic limit carefully. For this reason, we need to 
know asymptotic behavior of the graph of h{x) in the neighborhood of x = 0, which depends on the 
dimension of the system. 
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For D <2, h{0) tends to — oo as Q becomes large. Hence, for arbitrary (^^, x* is positive if Q is 
sufficiently large. Further, it is shown that 



X 



t > 



(29) 



because h{x) is strictly increasing. Combining eqs. (|25|) , (pBf) and (^Of) we conclude that uo{oo, ip^) is 
strictly convex and has the unique minimum at the origin. Therefore the order parameter vanishes 
in the thermodynamic limit. It is a typical feature of the effective potential in the symmetric phase. 
The result corresponds to the fact that a continuous symmetry is not broken when D <2 m & sense 
that a system has a non-vanishing order parameter. 

For -D > 3, /i(0) remains finite. If x > and sufficiently small, the integral in ( pH]) is evaluated 



as follows [23 



-C3 + TZ^ 



-C4 - 



47r 
1 



167r2 
'CD + KdX 



x\n.x 



{D = 3) 
iD = A) 
{D>4) 



where 



CD = 



{D>3), 



(30) 



(31) 



(27r)^-A(p)' 
and Kd is a non- universal positive constant. 

On the other hand, if —6 < x < 0, the graph of h{x) drags a tail as seen in Fig. 1. Since 5 ^ as 
Q — s> oo, the tail sticks along the y-axis, from (0, h{0)) towards (0, — oo) in the thermodynamic limit. 
This is the heart of the sticking argument pO|, |l^] used in the the spherical model [^. We emphasize 
that we could not have found the tail and have lost a real solution for + 2rn^/ < h{0) if we took 
the thermodynamic limit before solving the saddle-point equation |]T3|, p 



15, M 



From the above consideration for the asymptotic behavior of h{x), we conclude 



lim X* 



> {t'>^l) 
= it'<^l) 



(32) 



for D > 3. Here we defined the constant as 



2ml 
90 



- CD- 



(33) 



Note that ip^ can be both positive and negative. The first term of the right-hand-side in eq. (|3^) 
is related to the radius of the sphere in the target space where the classical potential takes the 
minimum. The second term originates from statistical fiuctuations. If v^o ^ 0' namely, if the 
statistical fiuctuations dominate over the classical radius, the sphere "disappears" and the 0{N) 
symmetry will be restored. 
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In fact, if < 0; then x* > for all t. Namely the effective potential in this case has the unique 
minimum at the origin and the symmetry is not spontaneously broken as in the lower dimensional 
case. 

In contrast, if ip^ > 0, x* vanishes when < ip^. This indicates that uo{oo, cp"^) is flat inside the 
(A^ — l)-dimensional sphere with the radius ipo in c/?-space according to eq.(^). This is an expected 
result in accord with numerical simulations |l]] or the Maxwell construction when the symmetry 
is spontaneously broken. The reason why the fiat region appears is that configurations with a domain 
structure contribute. [§, 0. 

We can quantitatively study the form of uo{oo, ip^) outside of the fiat region near the boundary, 
where we can assume (y?^ — ifl) is positive and sufficiently small. Since < x* << 1 in this case and 
the formulae (^) are available. For D = 3, we have 

X* = lQn\^'-^lf + -.., (34) 

where ■ ■ ■ stands for terms higher than (c/?^ — v^q)^- This provides, ignoring a constant term, 

I LZtt 
= + (35) 

Similarly, we have 

Using tio(oo,(^^), we can compute the one-point function in a usual manner]^, |T], |23|. First we 
introduce a uniform external field \fNJ while Vl is kept finite. That gives the following change in 
the partition function(^: 

[/o(^], (^2) ^ t/o(^^, ¥>') - Nr^\ (37) 

The one-point function is computed as 

(^E^) = lim I nrf^Vexp{-l] iV(«o(fi,c^') - JV^)}, (38) 

il ^ J^qU^oc Z[J\ J ^ 

where Z[J] is the partition function in the presence of the external field J. Note that the limit- 
ing procedure does not commute. Since the external field explicitly breaks the 0{N) symmetry, 
■uo(oo,<^^) — J°-(f"' takes the minimum at a unique point in c^-space, say, (p*. The integration over 
(f"' in (|38D can be exactly evaluated at ip* in the thermodynamic limit. As J 0, (p* approaches to 
(fQ e, where e is the unit vector pointing to the same direction as the external field. Thus we get 

{^j:^) = ^^oe. (39) 
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It shows that the symmetry is spontaneously broken. 

To summarize, in the case oi D > 3, the effective potential has a ffat region and the order 
parameter takes the value of (po e if and only if v^g > 0- Namely, the symmetry is spontaneously 
broken if and only if > 0. The direction e is determined by a uniform external field. 
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Fig. 2. The effective potential of a finite system in D = 4. 
Tfie classical potential(long-dashed line), = 2^ (dashed 
line) and f2 = 8** are presented. The coupling constants 
are chosen as go = 80, TOq — —15, which are identical with 
those of the Monte Carlo simulations in rcf. if we put 
iV = 1. 



3.2 Effective potential in a finite-size system 

In this subsection, we show numerical results for finite systems in order to see how the effective 
potential approaches to the form clarified in the previous subsection, as i7 — oo. We can numerically 
solve the saddle-point equation (|T3|) for finite systems and obtain Fig. 2. A similar result is obtained 
in the context of renormalization-group analysis p7|. We wish to understand, from a quantitative 
point of view, how the barrier of the effective potential flattens as the system size becomes large. To 
this end, we define the barrier height H as follows: 

if = Mo(0,l]) - minuo, (40) 

where we denote minwo as a minimum value of the effective potential. It is noted that an extreme 
value of the effective potential satisfies the condition = or x* = from eq. (^). When the 
system is in the broken phase, x* is negative at = for all fi, so that mq has a local maximum 
at 95 = 0. On the other hand, since h{x) of eq. (pTj) is strictly increasing as discussed above, there 
exists a solution x* = for some \ip\^ > 0, where uq has a local minimum. Thus, from eqs.(0) and 
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(|2^), the barrier height H is exphcitly written as 

1 ^ ^1 771^ 

H = 7^ T.M-m + X*m - ln(-A(p))} - —x*\0) + ^x*(0), (41) 

where 

X%0)^X*\ip=o. (42) 

Let us first consider how H depends on the coupling constants in a finite-size system. Since 
—6 < x*(0) < in the broken phase, x*(0) goes to zero as ^ oo. Therefore we can neglect the 
X*^(0)-term in eq. (^If) and also, due to the same reason, the x*-term in eq. (pISl). That means that 



H depends only on the combination m^/ Qq. Neglecting these terms can be effectively carried out by 
taking the limit qq ^ oo with (3 = —m^/gQ kept finite. In this limit the 0{N) 0^ model becomes the 
0{N) nonlinear a-model, where a length of a field variable is fixed to y/2N/3. 

We numerically compute the barrier height varying coupling constants with f3 kept fixed and 
plot the ratio of these barrier heights to that of the 0{N) nonlinear a-model as a function of L in 
D = 3(Fig. 3) and D = 4(Fig. 4). 



1.02 

0.98 ' = 
0.96 - 

20 40 60 80 100 ^ 

Fig. 3. The ratio of the barrier height H/Hq as a func- 
tion of L from L — 10 to 100 in £> = 3. The couphng 
constants are (ffoj^iQ) — (40, — 7.5)(boxes), (5oi"io) = 
(80, —15) (crosses), and (160, —30) (triangles), where P is 
kept at 0.1875. Hq is the barrier height of the 0{N) non- 
hnear cr-modeL 
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1.04 



1.02 
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20 
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Fig. 4. The ratio of the barrier height H/Hq as a function 
of L from L = 10 to 100 in D = 4. The couphng constants 
are the same as Fig. 3. 



The result shows that they have the same asymptotic behavior as that of the 0{N) nonhnear a-model 
for sufficiently large L, as we expected. Hence we conclude that, when the system is in the broken 
phase, the length of the field at each lattice is ^2N(3 in dominant configurations for sufficiently large 
L. 

In order to clarify dominant configuration at = more definitely, we treat if as a function of L 
and study its asymptotic behavior. When the parameters (3 and L are sufficiently large, contribution 
from energy to the effective potential will dominate over that from entropy in the broken phase. If 
the energy density (i.e., the value of the action per unit volume) of most dominant configurations at 
95 = behaves as 

min uq + const, x , (43) 

then H will indicate asymptotic behavior as const, x L~'^ . 

The barrier heights as a function of L in = 3 and D = A are shown in Fig. 5. We find that for 
L > 10 each H is fitted very well by the power function eq. (^3[) . The fitting exponent is summarized 
in Table 1. The result is that the all exponents are very close to 2. 
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Fig. 5. The barrier height H = uqCO, ^) — minwo a function of L from 
L = 2 to 100 with {go, ml) = (80,-15) in D = 3(daggers), {go, ml) = 
(80,-150) in L> = 3(boxes), {go, ml) = (80,-15) in L> = 4(crosses), 
{go, ml) = (80,-150) in D = 4(triangles). Each line is a fitting power 
function. 



D 


(f/o-"'o) 




3 


(80, -15) 


1.888 


3 


(80, -150) 


1.980 


4 


(80, -15) 


1.953 


4 


(80, -150) 


1.986 



Table 1: Fitting exponent 7 in the case of {go, ml) = (80, —15), (80, —150) in D = 3 and 4. 
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In order to understand the above asymptotic property, let us find a configuration satisfying the 
following conditions: 

1. c^ = 0, 

2. the length of the field at each lattice is y/2W^, 

3. its energy density behaves as const, x 

Let us observe a configuration where 0^ with length y/2Nj3 slowly rotates for one direction and 
is uniform for all the other directions in the lattice. For example^, 

0i = V^cos-^, 02^y2iV^sin-^, 0^ = (3<2<iV). (44) 

We find from eq. (|1]) that the energy density of this configuration measured from that of the uniform 
configuration = \/2Nf3 is A^/?(5, which asymptotically behaves as Hence such configurations 
satisfy the above conditions. 

It is noted that energy density of a configuration composed of some macroscopic uniform domains 
can be roughly estimated to be proportional to L^^. 

Therefore our numerical study strongly suggests that the barrier height of the effective potential in 
a finite-size system flattens as the system size becomes large because the slowly rotating configuration 
described above becomes dominant as the size becomes large. 

4 The continuum limit and renormalization effects 

So far, we have treated the effective potential of the lattice model. In this section, we take the 
continuum limit and study the renormalization effects. The global form of the effective potential of 
the continuum theory in dimension lower than three is obtained in refs.[^ |15|. Here we restrict the 
dimension of the system to three or four. Our main interest is to see a fate of the fiat region. 

In order to take the continuum limit, we introduce the lattice spacing parameter a, which measures 
dimensions of physical quantities. We define 

V = a^VL, 

b -2^ b 

I, D-2 L 

0^ ^ a- — 0^ 

Xr = (45) 

•^Similar configuration was considered by Ringwald and Wetterich[p|. 
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where we assign the canonical dimension to the scalar fields. The effective potential of the lattice 
theory with the finite volume V and with the lattice parameter a, u^{a, V, 95^), is defined by 



UKia,V,cpl) = a-''uoin,ip-'). (46) 



It is explicitly written in terms of dimensionful quantities defined in eq.(^). Using eqs. ([T^) and (|2 

we get 

unia, V,vl) = Iflxl - ^^^(^^X? + J^Xr + ^ E lu(-A(pa)a- + xkl (47) 

Here Xr is the solution of the saddle-point equation of the theory with the lattice spacing a, which 
is derived from eq.(|TE|): 

where the momentum is quantized by and runs from — vr/a to vr/a. We adopt a renor- 

malization prescription to keep MR(a, V, 95^) finite in the continuum limit a ^ 0. This determines 
lattice-spacing dependence of the bare parameters, gQ^a) and mQ^a). 

We shall first consider the case of D = 3. The bare parameters are chosen in such a way that the 
following equations hold|]T3|: 

9k = go{a)a~^, 
2m| ^ 2mg(a)a-^ 1 ^ 1 
gn go{a)a-^ V f^^ -K{pa)a-^' 

where we introduced the renormalized mass and the coupling constant, m\ and g-^. They are directly 
related to physical observables of the continuum theory. The continuum limit a — is taken keeping 
V^R5 ^Jii 9^.1 and V finite. In this limit, we can make the following substitution in (|47|) and (^) 

- A(pa)a-2 + Xk^P'' + Xr, (50) 

where = J2^PfiPtJ.- The continuum effective potential is written in terms of the renormalized 
coupling constants 



«r(0, V, cpl) = IcplxK - ^Xr' + ^Xr + ^ T.Mp' + Xr) - -,Xk}, (51) 
2 ^gn gn -^^ p 



where Xr satisfies 



2ml I 1 1 . 2 
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MR(00,t2) 




Fig. 6. The effective potential of the continuum theory in 
D ~ 3 and the classical potential. The renornialized coupling 
constants are chosen as — —15 and = 80. 



Our aim is to clarify the global form of the effective potential in the infinite- volume limit. Behavior 
of the solution in the limit V — oo is the same as that in the lattice model. Namely, if > 0, 



hm Xk > 0, (53) 

V—i-oo 



lim Xki" " ""1'''"'' . (54) 



and if < 0, 



We conclude from this behavior that mr(0, oo, </Jr) is flat inside the sphere with the radius \J —2m\/ g-p, 
when the symmetry is broken {m^ < 0), as is the case of the lattice model. When x*r > 0) the 
summation over p in (|52|) is replaced by the integration and the result is 

5 2m?. I I — 2 , 

We can explicitly solve the equation. Inserting the solution into (0), we obtain mr(0, oo, 9?r) outside 
the fiat region. We depict the result in Fig. 6. 

It may possible to take other limiting procedures for obtaining mr(0, oo, ^p^. For example, one 
may first take the limit Vt ^ oo with fixed a > and next take the limit a ^ 0. It corresponds to 
taking the infinite-volume limit before taking the continuum limit. We can obtain the same result in 
this limiting procedure if we solve the saddle-point equation before taking the limit Vt ^ oo. 

Let us turn to the case oi D = A. We introduce the renormalized coupling constants by the 



following equations [|13 



2 _ 2 1 ^ 1 
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2m| ^ 2ml{a)a ^ 1 1 /^gx 

where is an arbitrary positive constant. The saddle-point equation can be written in the finite 
form as 

2 2m| 1 1 12^ 

^""^l^^ V^,f^WT^ " WTW^^ + ^^^^ 



Since the summation over p in the first equation in eq. ([5B|) diverges as a — > 0, must be chosen 
to zero as long as we start from a positive (7o(ct)[|ll) IHl- This result corresponds to the conjecture 
that the (pf theory becomes trivial in the continuum limit. The saddle-point equation becomes 

< = Xk- (58) 
Namely, the theory reduces to the Gaussian model. The effective potential in the continuum limit is 

ur{0, V, ifl) = Iml^l (59) 

up to an irrelevant additive constant. We must choose m'^ to be non- negative for well-defined 
continuum filed theory. 



5 Summary and discussion 

In this paper, we obtained the effective potential of the lattice 0^ theory in the large- limit. In the 
thermodynamic limit, we found out that it is globally real- valued and convex in both the symmetric 
and the broken phases even if classical potential is not convex. In particular, the effective potential 
has a flat region in the broken phase, which is consistent with numerical simulations and with the 
Maxwell construction. 

The effective potential in a flnite system has a potential barrier at the origin although it is much 
lower than that of the classical potential. The height of the barrier is related to energy of defects 
of non- uniform conflgurations, and asymptotically flattens as L~^, where L is the linear size of the 
system. The result concerning with a flnite lattice system can be applied to checking the validity of 



an expansion p9[|. The effective potential with sufficiently flat maximum as shown in Fig. 2 may 



shed hght on global structure of inflationary universe pO[ . 

We also studied the continuum limit of this model and discussed effects of the renormalization 
to the flat region in the case of D = 3 and Z) = 4. In the case of D = 3, the effective potential in 
the continuum limit has the flat region as in the case of the lattice model. On the other hand, in the 
case of D = 4, the flat region disappears because the (j)'^ couphng is renormalized to zero. 
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Here we discuss a role of the effective potential in continuum quantum field theory. After taking 
the continuum limit, it is believed that a minimum of the effective potential corresponds to the 
ground state (vacuum) in quantum field theory]^. If there are multiple points giving the minimum, 
we need to choose one point. A standard way of doing this is to turn on a uniform external field 
as we demonstrated in sect. 3. The vacuum selected in this way is homogeneous and corresponds 
to a point on the boundary of the flat region. Since it is the vacuum that is used in the previous 



literature |jT3|, |T4|, |T6|, we can rederive the one-point function obtained in those works. However, 
physical observables concerned with higher-point correlation functions such as a mass matrix (or a 
susceptibility tensor) may be changed if we take the flat region into account. We hope to report on 
this point more detail in a subsequent publication. 

In addition to reproducing the previous results, the effective potential obtained here suggests 
possibility of a non-uniform ground state such as a rotating configuration with a long-wave length 
in eq.(l^) or as kink ground states in the quantum ferromagnetic XXZ chain[Q. In order to study 
such a non-homogeneous vacuum, it is expected that a non-uniform external field will be needed, 
which will pick out a point inside the fiat region. 

Related with this topic, we comment on the expansion in the presence of a non-uniform 
external field Jx- In this case, the action of the system changes to S" — \/NY,x Jx't'x- The strategy 
adopted in sect. 2 leads us to the following saddle-point equation 

1 1 1 

- Y.iiJ'^ + A-\n^'^ - i?")}{ J: + A-\n^'^ -B'^)}- A-')Gx-.G,y + -G,, - -X, + ^ = 0, (60) 

^ x,y ^ 90 90 

where 

Gxy^{x\—^\y), A^Y^Gxy, B'^^Y.Gxy.ry. (61) 

x,y 

Although we need to solve this equation for carrying out the 1/A^ expansion, it is difficult because 
we cannot employ the ansatz of a uniform solution due to the first term Q In fact, using the ansatz, 
we find that the constant solution x* must satisfy 

^ E G.{p) - -X* + ^ + (^ E i"(p)Go(p)e-)^ = 0, (62) 

for an arbitrary site z. Here S"'{p), P 7^ are identical with Fourier components of and 

J"(0) = fixV"- (63) 
However it is obviously impossible if an oscillating mode of Jx is present. 



ref. where a non-uniform external field is employed as a "regularization" , the first term is oversimplified and 
it misleads to a uniform solution. 
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Moreover, a perturbative expansion in the external source fails even if the magnitude of it is 
infinitesimal because the expansion induces the infrared divergence in the broken phase. 

We feel that a constant mode of (px is not appropriate for studying a non-uniform ground state. 
The order parameters cannot resolve whether an inner point of the flat region corresponds to an 
inhomogeneous states or a superposition of pure states selected by a uniform external source. Some 
oscillating modes will be chosen as order parameters. 
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